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Abstract. We report on a near- infrared adaptive optics 
survey of a sample of 66 low-mass members of the pre- 
main sequence stellar cluster IC348. We find 12 binary 
systems in the separation range 0'.'l-8", excluding 3 proba- 
ble background projected companions. An estimate of the 
number of faint undetected companions is derived, before 
we evaluate the binary frequency in this cluster. In the 
range logP = 5.0-7.9 days, the binary fraction in IC348 
is 19 ± 5 %. This is similar to the values corresponding to 
G- and M-dwarfs in the solar neighbourhood population 
(23 ±3 % and ~ 18%, respectively). Furthermore, the dis- 
tribution of orbital periods of IC 348 binaries in this range 
is consistent with that of field binaries. We conclude that 
there is no binary excess in IC 348. 

Substellar companions are found to be rare, or even 
missing, as companions of low-mass stars in the separation 
range we surveyed. Also, the mass ratio distribution is 
not peaked at q ~ 1 in IC348, and it is unlikely that an 
observational bias can account for that. 

We do not find any evidence for an evolution of the 
binary frequency with age within the age spread of the 
cluster of about lOMyr. Comparing the binary frequency 
in IC 348 with that of other star forming regions (SFRs) 
and young open clusters, we conclude that there is no sig- 
nificant temporal evolution of the binary fraction between 
a few Myrs after the formation process and the zero-age 
main sequence (ZAMS) and field populations. We find in- 
stead a trend for the binary fraction to be inversely cor- 
related with stellar density, with dense clusters having a 
binary fraction similar to that of field dwarfs and loose 
associations exhibiting an excess of binaries. Two scenar- 
ios can be suggested to explain these differences: either 
all SFRs, clusters and associations alike, initially host a 
large number of binaries, which is subsequently reduced 
only in dense clusters on a timescale of less than 1 Myr 
due to numerous gravitational encounters, or specific ini- 
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tial conditions in the parental molecular clouds impact on 
the fragmentation process leading to intrinsically different 
binary fractions from one SFR to the other. 

Key words: binaries: visual - stars: pre-main sequence - 
stars: imaging - stars: late-type - stars: formation 



1. Introduction 

Several studies in the early 90s have shown that binarity 
is a very common property of low-mass main sequence 
stars: about 53% of G-type stars, 45% of K-dwarfs and 
42% of M-dwarfs are in fact multiple in the solar neigh- 
bourhood (Duquennoy & Mayor 1991, hereafter DM91; 
Mayor et al. [Tgggt Fischer & Marcy |1992|). An important 



issue for current star formation models is to account for 
the high number of binaries, and to predict their physical 
properties. 

The relative number of binary systems may be even 
larger among much younger stellar populations. One of 
the best studied low-mass SFRs, the Taurus- Auriga dark 
cloud, hosts almost twice as many binaries as the solar 
neighbourhood in the separation range 2-2000 AU (Lein- 
ert et al. |1993[ Ghez ct al. |1993| ; Simon et al. |1995[ ). Yet, 
subsequent surveys of a number of other SFRs have led 
to somewhat conflicting results: while some exhibit binary 
excesses comparable to that of the Taurus cloud (Padgett 
et al. 1997; Ghez et al. 1997), others appear to have simi- 



lar binary fractions as field dwarfs (Brandner et al. 1996 



For instance, the Orion Trapezium shows a binary fraction 



in good agreement with that on the MS (Petr et al. 1998 



Prosser et al. 1994). Duchene (1998) recently reanalysed 



in a consistent way these various studies and confirmed 
that the binary fraction appears to vary from one SFR to 
the other, with the main exception of all Orion clusters 
whose binary fraction is similar to that of the field. 

Several proposals have been made to account for these 
results. It has been suggested that the fragmentation pro- 
cess during the protostellar collapse yields a high fraction 
of multiple systems which, however, steadily declines over 
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time as multiple systems are disrupted during their sub- 
sequent evolution (Ghez et al. 1993). Then, the binary 



fraction would depend on the age of the stars and would 
vary over a timescale of several 100 Myr (Patience et al. 



1998). Alternatively, it is conceivable that the binary frac- 



tion of a cloud is established at the very beginning of the 



cluster history. Kroupa (1995a, 1995b) has recently shown 
from TV-body simulations that in regions as dense as the 
Trapezium cluster the binary fraction could decrease from 
100% to about 50% in less than 1 Myr due to gravita- 
tional encounters. Still another possibility is that the bi- 
nary frequency is sensitive to environmental conditions 
in the parental molecular cloud. In a qualitative study, 
Durisen & Sterzick ( |l994|) found that both the current 
fragmentation models and disk instabilities are compat- 
ible with a lower binary fraction in clouds with higher 
temperatures. 

In order to distinguish between these alternatives, we 
started a long-term project aimed at studying binaries in 
clusters at different evolutionary stages, from the birthline 
to the MS. In Bouvier et al. ( 1997 ), we already found that 
the binary frequency in the 100 Myr old Pleiades cluster 
is similar to that of the MS. We report here the results 
obtained in the pre-main sequence (PMS) cluster IC348. 
This cluster was selected on the basis of its age being 
similar to that of the Taurus cloud (about 2 Myr), but 
its stellar density being much larger (about 500 stars pc -3 
compared to a few stars pc -3 for Taurus). 

IC 348 is a young cluster located in the Perseus molec- 



ular cloud, at a distance of about 320 pc (Herbig 1998, 
hereafter H98). It hosts a B5 V star (BD +31°643), about 
100 optical sources (Trullols & Jordi 1997), as well as a few 



hundred infrared sources, which are probably embedded 
young stars (Lada & Lada 1995 ). The age of this cluster 
has been estimated by several methods: Lada & Lada, by 
fitting the IMF, estimate that star formation is s till go - 
ing on after a burst ~ 5-7 Myr ago; Luhman et al. ( 1998[ ), 
hereafter L98, found that a major burst occured ~ 3 Myr 
ago, but that stars as old as 10 Myr also lie in the cluster. 
Similarly, from a dereddened colour-magnitude diagram, 
H98 estimates ages ranging from less than 1 Myr and up 
to about 10 Myr for about 100 members, with a median 
age of ~ 2 Myr. From a deep near-infrared survey, Lada & 
Lada estimated a stellar density of about 500 stars pc -3 or 
220M Q pc -3 (within the half-mass radius of 0.47 pc), and 
a projected surface density of about 1000 stars pc -2 in the 
central 0.1 pc, similar to the NGC 2024 cluster in Orion. 
H98 also evaluated masses for the members, and found 
that the median mass is about 0.5Mq, in agreement with 
the IMF estimated for this cluster (which is very similar to 
the IMF from Scalo 1986). Furthermore, H98 conducted 
an Ha survey, and discovered over 110 emission- line stars 
in a 7.'5 x 15' area centred on the cluster; all of these stars 
are very likely to be young, active cluster members; about 
70 of them have independently been confirmed as mem- 
bers on the basis of colour-magnitude diagrams. Preibisch 



et al. (1996) performed a systematic X-ray survey of the 
area, detecting over 110 sources within a l°-radius circle. 

We report on our adaptive optics observations of 
IC 348 in Sect.B, and estimate the binary fraction in the 
cluster in Sect. [| In Sect. || we discuss various other binary 
properties (especially the scarcity of brown dwarf compan- 
ions in our survey), and Sect.^j presents a discussion on 
the link between the binary fraction and environmental 
conditions in SFRs. Sect.|^ summarizes our conclusions. 

2. Observations 

Our sample was selected from H98's list: all stars brighter 
than R « 15 were observed, with a 75% completeness 
level at R = 16.5; overall, the survey is two-third com- 
plete at the H = 12 limit. Mosaicing with small offsets 
from the brightest stars, about 25 additional fainter stars 
were surveyed. In total, we observed 70 of H98's mem- 
bers (within 66 independent systems); we also observed 
24 stars for which H98 could not assess membership. Ten 
of these stars are considered as members by L98, on the 
basis of the detection of the Li A6707A absorption line 
and of spectral classification (for late M stars). 

The observations were obtained during four nights in 
December 1997 at the Canada-France-Hawaii Telescope 
on Mauna Kea. We performed near-infrared (JHK), high 
angular resolution imaging using the Adaptive Optics 
Bonnette and the new infrared camera KIR, a 1024 x 1024 
HAWAII detector. The pixel scale is 0'.'0351/pix, yield- 
ing a total field-of-view of 36" x 36". Most of the images 
are diffraction-limited in H and K, but the images of the 
faintest stars or those observed at large offsets from the 
wavefront star (> 20"), have FWHM as large as 0"2; Fig.[| 
illustrates the image quality in our survey. We surveyed 
each target in H, and all resolved systems were also ob- 
served in J and K (with the exception of pairs formed 
by two stars already known as members). A typical ob- 
serving sequence consists of 16 images, at 4 positions, 
with individual exposure times ranging from 2 to 30 sec- 
onds to avoid saturation of the sources. On the first night, 
the observing procedure was different, since we produced 
1 '5 x 1 '5 mosaics centred on the brightest stars close to 
the cluster centre. In these images, the exposure time was 
such that no star was saturated; the 5 a detection limit 
in these images is about H = 15.5. For all but two H98 
members (If A 134 and 163), this leads to a detection limit 
of at least sa 3 mag for separations larger than 0"5. 

UKIRT photometric standards were observed every 
night for the flux calibration. Data reduction was per- 
formed with IRAF packages. Astrometry and relative pho- 
tometry of binaries were obtained by PSF fitting and then 
combined with large-aperture photometry to get the abso- 
lute photometry of each component. In two cases (If A 100, 
102) deconvolution had to be applied to obtain the rel- 
ative photometry; we used the Iterative Deconvolution 
Analysis in C (IDAC) routine (Jefferies & Christou 1993 ), 
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(a) 



(b) 



(c) 



(d) 



Fig. 1. Examples of images obtained in our survey. Each box is 1'.'75 x 1'.'75. (a): IfA 166 is a bright star (R — 11.2), 
allowing diffraction-limited images, with high Strehl ratios at K - (b): the fainter (R = 15.8) binary IfA 119 is still 
resolved at K with a separation of 0'.'25 - (c) IfA 192 is too faint (R — 16.5) to be resolved as a binary (separation 
of 0'.'13), but it is clearly elongated at J, H and K, and a deconvolution process was used to obtain the relative 
photometry - (d): IfA 184 illustrates the limit for the detection of close, faint companions (the image is in the J band). 



, back. 



separation (") 



Fig. 2. The observed binaries in our sample, including the 
probable background stars (labelled "back."). Open circles 
represent L98's members detected as binaries. The error 
bars are smaller than the symbol size, except for IfA 184 
(see text). The dashed line is an estimate of the detection 
limit: it is the 3cr noise level in the averaged radial pro- 
file of a single star. The histogram is another estimate of 
this limit, obtained by artificially adding faint companions 
close to single stars. 



and cross-checked with the Lucy algorithm. In the case of 
IfA 184, deconvolution was unsuccessful in the H band, 
and we consider here the results from PSF fitting, al- 
though we think that this method leads to an underes- 
timate of the flux ratio in this specific case. Estimates of 
the uncertainties are 0.05 mag and 0.02 mag for absolute 
and relative photometry, 0'.'005 for the separation and 0?2 
for the position angle. The errors are slightly larger when 
deconvolution was applied; when the adaptive optics sys- 



tem was locked on a binary system (e.g., IfA 139-140), 
the PSF was substantially deformed, leading to increased 
uncertainties on the centroid locations. 

3. Binary fraction in IC 348 

Our binary candidates are listed in Table[l], together with 
their astrometric and photometric measurements. To ac- 
cept two stars as a binary candidate, we have set an up- 
per limit for the separation of 8", corresponding to ^2500 
AU. This limit should avoid confusion between real com- 
panions and background stars (see discussion below). A 
case-by-case study has been done, however, for each bi- 
nary candidate. All the stars that appeared as singles are 
listed in Table|. The astrometry for the binaries detected 
by H98 in the / band agrees with ours to within 0'.'15 and 
2°, even for the closest pairs (IfA 136 and 211). We failed 
to detect the faint companion 1'.'3 away from IfA 159, be- 
cause the latter is a faint star (R = 18.7, H = 12.19), 
which was observed in a mosaic on the first night, with 
a more limited dynamic range compared to that usually 
achieved on other stars. At the 5 a level, the companion is 
fainter than J — 16.2 and H — 15.2 mag or, equivalently, 
the binary flux ratio is larger than 2.9 mag in both bands. 

Fig. [2] shows the magnitude difference in the H band 
as a function of separation for detected binaries. The solid 
and dotted lines show the detection limit of our survey, 
estimated in two different ways: the solid histogram was 
established by adding faint stars around single targets 
and visually inspecting the images. Since it corresponds 
roughly to a 5 a peak detection, it lies 0.5 to 1 mag above 
the dotted line, which represents the 3 a noise level as mea- 
sured on the PSFs of single stars in our images. At large 
distances from the primaries (> 1"), companions can be 
detected down to ~ 6.5 mag fainter than the primaries. In 
some cases, even somewhat fainter stars can be detected. 
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Table 1. Astrometry and photometry for binary candidates in IC348. Position angles are given East from 
North. Boldface entries in the first column are H98's members. The primary mass is from Herbig (priv. comm.), 
and takes into account individual ages; the mass-ratio is estimated from AH and the mass-luminosity relation 
M H = -3.251og(M/M©) + 2.19 for M < 1M , Baraffe et al. ( |1998| ). Both methods are somewhat uncertai n but 
are independent. TJ 81 =BD+31°643 does not have an IfA number, but is designed by its Trullos & Jordi (1997) 
identification. 
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M B /M a 


Probable member binaries 


48-49 




12.41 






0.99 




1.343 


112.4 


0.2 


0.50 


85-82 




11.70 






1.89 




3.086 


289.2 


0.4 


0.26 


102 1 


12.29 


11.25 


11.10 


0.69±0.05 


0.67±0.05 


0.62±0.05 


0.26±0.01 


357.0±0.5 


0.4 


0.62 


119 


12.16 


11.30 


11.00 


0.63 


0.67 


0.65 


0.249 


254.0 


0.3 


0.62 


136 2 


13.43 


12.63 


12.33 


0.74±0.04 


0.75±0.04 


0.71±0.04 


0.70±0.03 


214±1 


0.2 


0.59 


139-140 1 " 2 




9.17 






0.26±0.04 




1.23±0.03 


86±1 


1.4 


0.83 


144-143 


9.99 


9.21 




1.01 


0.73 




5.371 


338.3 


1.4 


0.60 


157-158 


10.32 


9.52 


9.01 


0.95 


0.72 


0.46 


4.321 


101.0 


0.5 


0.60 


166"^ 


9.09 


8.33 


8.11 


2.04 


1.49 


1.24 


0.559 


151.3 




0.35 


184 3 


9.95 


9.07 


8.70 


2.55 


l-84_ 10 


2.08 


0.409 


349.7 


0.4 


0.27 


192 1 




12.97 


12.41 




1.35±0.05 


0.87±0.05 


0.13±0.01 


186.6±0.5 


0.2 


0.38 


211 


12.23 


10.84 


10.46 


1.52 


1.48 


1.40 


1.028 


155.6 


0.7 


0.35 


261-104 




8.63 






1.46 




6.529 


312.6 


1.9 


0.36 


TJ81 n 


6.76 


6.53 


6.51 


0.20±0.04 


0.14±0.04 


0.14±0.04 


0.61±0.03 


16±1 




0.91 


Probable background companions and non-member primaries 


20-21 




9.63 






1.88 




2.129 


112.4 






100 1 


14.63 


13.71 


13.49 






0.62±0.05 


0.13±0.01 


100.7±0.5 






124 4 




10.92 






5.38±0.07 




2.536 


103.2 






137 4 


10.63 


9.86 




5.86±0.09 


5.03±0.08 




6.734 


314.5 






LkHa86 4 




11.03 






7.44±0.09 




2.588 


339.2 







1 deconvolved images 
as members by L98. 



AO locked on a binary system - 3 poor PSF fitting in H - 4 background companion - ' stars identified 



Table 2. Stars that are unresolved in our survey; boldface entries are for H98's members. 
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stars identified as members by L98. 



The fainter stars observed in the mosaics on the first night 
were observed with smaller signal-to-noise ratios. The de- 
tection limit at large separations is thus poorer for these 
stars than in Fig-H; close to the primaries, however, the 
detection limit remains roughly unchanged, since the lim- 
itation comes primarily from photon noise. 

Despite a large dynamic range in our images, we found 
only 3 secondaries with AH > 2.5, with in fact AH > 5. 
The location of the widest of these companions (marked 
"back." in Fig.|J) in a J-(J — H) diagram indicates that 
it is a background star, lying well away from all known 
cluster members. The two other very faint companions 
(marked "back.?") are also likely background stars, al- 
though we lack multicolour photometry to prove it. For 
the binary IfA 139-IfA 140, we only have H photometry 
available; we believe it is a physical binary, however, be- 



cause of its close separation. The membership of IfA 140 
could not be determined by H98 since no / band pho- 
tometry was obtained for this star; L98 classify it as a 
member. Similarly, IfA 82 lacks V measurement in H98's 
study, so that its membership is not decided. We consider, 
however, the pair IfA85-IfA82 as a physical binary, be- 
cause of the late spectral type of IfA 82 (M4, H98). For 
all other systems, the location of the companions in a H— 
(H — K) diagram suggests membership and thus physical 
association with the primary. 

Close to the primaries, it is difficult to detect very faint 
companions, because of the wings of the PSF: below 0'.'2, 
only binaries with flux ratios AH < 2 can be resolved. 
To calculate the actual binary fraction in IC 348, we need 
to estimate the number of fainter, undetected secondaries 
at these separations. The method we use, fully described 
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Table 3. Completeness correction for this survey, "fraction missed" is the ratio of the number of undetected companions 
to the total number of companions. The last column summarizes the overall figures. Only H98 members are considered 
here. 



sep. range (") 


0.1-0.2 


0.2-0.3 


0.3-0.5 


0.5-0.75 


0.75-1.0 


1.0-2.0 


2.0-8.0 


0.1-8.0 


orbital period (log(d)) 


5.0-5.4 


5.4-5.7 


5.7-6.0 


6.0-6.3 


6.3-6.5 


6.5-6.9 


6.9-7.8 


5.0-7.8 


AH Um (mag) 


2.0 


2.5 


3.0 


4.0 


5.0 


6.0 


6.5 




Qmin 


0.24 


0.17 


0.12 


< 0.1 


< 0.1 


< 0.1 


< 0.1 




fraction missed 


22% 


9% 


3% 










4% 


detected companions 


1 


2 


1 


1 





3 


4 


12 


corrected companions 


1.3 


2.2 


1.0 


1 





3 


4 


12.5 


IC 348 binary fraction (%) 




6.8±3.4 






6.1±3.0 




6.1±3.0 


18.9±5.3 


G-dwarf fraction (%) 




10.6±1.2 






7.4±0.8 




5.2±0.6 


23.1±2.6 



in Bouvier et al. ( 1997 ) for their survey of the Pleiades, 
consists in estimating the detection limit in several sepa- 
ration bins (chosen such that it is roughly constant within 
each bin); this flux ratio limit is then converted into a 
limiting mass ratio using the mass-luminosity relationship 
for 2Myr-old stars from Baraffe et al. ( |1998| ). Finally it 
is assumed that the mass ratio distribution observed by 
DM91 in the solar neighbourhood for G-dwarf binaries 
applies to IC 348 binaries. Then, the limiting mass ratio 
can be transformed into a fraction of missed companions. 
Because of the young age of the cluster: a mass ratio of 
q = 0.1, which is the limit of the DM91 survey, corre- 
sponds to AH ~ 3 mag at 2 Myr. It can be seen, from 
Fig.Q, that we reached this flux ratio over almost the en- 
tire range of separation considered. For the innermost Of.' 3, 
where the limiting flux ratio of our observations is slightly 
smaller, companions close to the q = 0.1 limit remain un- 
detected. Therefore, the overall correction is small, with 
only about 4% of the companions missed. We estimate 
the number of such companions per bin of separation in 
Table |. The mass ratio distribution of M-dwa rfs m ay be 
flatter t han tha t of G-dwarfs (Fischer & Marcy 1992 ; Rcid 
& Gizis 1997a ), leading to a different estimation of the 
completeness correction. However, a flatter distribution 
implies that we have missed even less companions (the 
number of binaries with small q is smaller), so that our 
estimate of the number of missed companions can be con- 
sidered as a conservative upper value. 

In Table |[ we also estimate the MS binary fraction 
(" G-dwarf fraction" ) in different separation ranges by in- 
tegrating the binary distribution from DM91. This is the 
number of binaries as a function of the orbital period. 
For IC348 binaries, we only have the angular separation 
between both components. Therefore, we need to convert 
these separations into orbital periods. We use the distance 
to the cluster and a statistical correction for the projec- 
tion of the semi-major axis on the sky (log a = log/? + 0.1, 
where a is the actual semi-major axis and p is the ap- 



parent separation, Reipurth & Zinnecker 1993). We also 
assume that the average total mass of a system is 1M Q 
(mean value for the observed binaries); a typical mass of 



5M© was assumed for BD+31°643. Small changes in the 
assumptions about the distance to the cluster, the projec- 
tion correction, and the stellar masses do not change the 
results significantly, since the MS orbital period distribu- 
tion is very broad: a small shift in the integration bound- 
aries does not significantly modify the binary fraction in 
this range. 

The overall binary fraction (number of companions per 
observed target) in IC 348 is 19 ± 5 % in the separation 
range 40-3200 AU; here, linear separations have been cor- 
rected to account for projection effects. It is not modified 
if L98's members are included. From DM91, we evaluate 
that the MS G-dwarf binary fraction over the same range 
is 23 ± 3 %. M-dwarfs have a binary fraction of about 18% 
over the same separation range (Fishcer & Marcy 1992), 
again very similar to IC 348. We checked that the presence 
of a few stars with early spectral type (earlier than K0) in 
our sample does not bias our results: the binary fraction 
in the subsample of stars with spectral type later than K0 
is 17 ± 7 %, indistinguishable of that of the whole sample. 

The orbital period distribution of IC 348 binaries is 
shown in Fig. |[ The comparison with the MS for each bin 
of the histogram is also given in Table ||. Both the plot 
in Fig.|3| and the similar values for the binary fraction in 
the three separation ranges in Table || demonstrate that 
the observed distribution is rather flat and, within the 
errors, not different from the MS distribution. It is again 
noticeable that the results are not strongly modified if we 
include or exclude higher mass stars or L98's members. 
We will now only consider H98 members, for which masses 
and ages have been determined. 

To evaluate the impact of possible background stars, 
we estimated the binary fraction over a smaller separa- 
tion range, with an upper limit set at 2" (640 AU). In 
this case, the binary fraction for IC 348 and the MS are 
respectively 13 ± 5 % and 19 ± 2 %, i.e. similar values, 
possibly with a small deficiency in IC348. This suggests 
that the number of false detections is small in the sepa- 
ration range we selected. On the other hand, in the range 
8-16", we find 6 "companions", that is a binary fraction 
of 9 ±4 %, while the MS value is 1.8 ±0.2 %. Furthermore, 
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any triple or higher order multiple system. Then, the ac- 
tivity properties of binaries are compared to those of single 



log P (days) 



Fig. 3. Orbital period distribution in IC 348 (solid his- 
togram), compared to the empirical distributions for G- 
(dotted c urve, DM91) and M-type (dashed line, Fischer 
& Marcy 1992 ) MS stars. The long-dashed histogram in- 
cludes L98's members, while the dotted-dashed histogram 
represent the subsample of stars with spectral type later 
than KO. For clarity, error bars are only drawn for one 
histogram. 

the only companion in this separation range for which we 
have two-colour photometry appears to be a background 
star. This supports our choice for the upper limit: larger 
values would imply a non-negligible background star con- 
tamination (unless the period distribution were different 
from that in the MS, with a peak at much larger separa- 
tions). Using the st ellar d ensity in the H off- field image 
from Lada & Lada (1995) without any correction for ex- 
tinction (i.e., overestimating the number of background 
companions), we expect about 4 false detections, similar 
to our findings. On the other hand, one has to determine 
the occurence of projection pairing of two members. The 
projected stellar density of the cluster can be crudely esti- 
mated from Lada & Lada ( 1995| )'s survey. Once field star 
contamination is substracted, they are left with an aver- 
age density of 2.5 10 -4 stars/" 2 , implying a total of about 
3 members chance projection in the range Qf.'l— 8". The 
pair IfA261-IfA 104 is a good candidate for such a pro- 
jection effect, since two other members (IfA 106 and 107) 
lie within 15" away from IfA 261. We do not try to correct 
for this effect since we lack a local estimate of the stellar 
density around each binary candidate. 



4. Binary properties in IC 348 

Various properties of stars, especially of PMS objects, de- 
pend on the stellar environment. The presence of a com- 
panion in the vicinity of a star modifies this environment 
in a non-negligible way. Potentially, this can affect the 
physi cal properties of the stars in multiple systems. In 
Sect. 4.1 we first consider the fact that we have not found 



stars (Sect. 4.2). Finally, our non-detection of very small 
mass ratios (q < 0.25) is discussed in Sect. 4.2, as well as 
the absence of candidate brown dwarfs. 

4-1. High order multiple systems 

In our survey, we have found 12 binaries, but no triple or 
quadruple systems. Does this mean that there is a defi- 
ciency of higher order multiple systems in IC 348? 

In the solar neighbourhood, G-dwarfs host roughly 10 
binaries for 1 triple system, and 4 triple systems for each 
quadruple system (DM91). Systems with more than two 
stars are thus quite rare. For Taurus PMS stars, Leinert 



et al. (1993) find a ratio of binaries to higher order mul- 
tiples of about 8:1, showing that the number of triples 
and quadruples does not seem to evolve significantly from 
PMS to MS stages. If we assume that this ratio of 10 to 1 
is also relevant for IC 348, then we would have expected to 
find one triple system, which is not statistically different 
from our findings, given the small numbers involved. 

There may be a second reason why we failed to detect 
triple systems. These systems are usually hierarchical (in 
both MS and PMS populations), with a close system sur- 
rounded by the orbit of a third star lying further away. 
Usually, the ratio of the two semi-major axes in triple sys- 
tems is at least «5 (Tokovinin 1997). Given the distance 



to IC 348, the peak of the orbital period distribution cor- 
responds roughly to the smallest separation we can resolve 
(see Fig.||), and most of the triple systems that have been 
detected in other SFRs have their orbits on both sides of 
this peak. Therefore, in IC348, triples may just appear 
as wide binaries, with the close binary system remaining 
undetected. 

The absence of any triple system from our sample most 
likely is not an indication for a different binary-to-triplc 
ratio between IC348 and the MS; it is probably due to 
the distance to the cluster and to statistical uncertainties. 

4-2. Stellar activity in binary systems 

Several indicators of T Tauri star activity have been iden- 
tified so far. X-ray and Balmer line (e.g., Ha) emission are 
some of these indicators. We have searched in our sample 
for a possible impact of binarity upon this activity (see 
Table ||). Although the samples are rather small, it ap- 
pears that binaries and single stars have roughly the same 
fraction of emitting stars. 

It seems, from our survey, that the magnetic activ- 
ity of PMS stars is not dependent on the presence of a 
companion. Binary surveys in X-ray selected samples had 
to face the bias induced in their target selection by the 
fact that both stars can be emitters. Various estimates 



of this bias were o btaine d by Brandner et al. ( 1996 ) and 
Kohler & Leinert ( 1998 ), but it appeared that it is not 
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Table 4. Comparison of the overall and binary samples 
regarding their X-ray and Ha emission. A binary is con- 
sidered as emitting if at least one of its components shows 
emission. 







observed 


binaries 


ROSAT 


sources 


30 


6 




undetected 


36 


6 


Ha 


emission 


31 


6 




unknown or abs. 


35 


6 



an important effect. This is confirmed by the absence of 
a significant difference in the binary fraction of ROSAT- 
detected and undetected stars in a single cluster (Leinert 
et al. 11993} Kohler & Leinert |1998|). Similarly, in IC348, 



the X-ray source sample does not show a higher binary 
frequency than the whole sample. 

Because Ha emission is likely linked to the accretion 
phenomenon on TTauri stars (e.g., Edwards et al. 1994), 
it also appears that binary members in IC 348 are sur- 
rounded by accretion disks in the same proportion as sin- 
gle stars, i.e., that binaries do not disrupt disks more 
rapidly than singles. Indeed, if two stars are separated 
by a few tens of AU, inner disks can remain unaffected 
around these stars. This is enough to support accretion 
onto the stars and to emit Balmer lines. Similarly, Prato 



& Simon (1997) showed that the near- infrared emission 
of accretion disks is the same in multiple systems and in 
single TTauri stars in the Taurus- Auriga SFR, provided 
the companions are separated by 40 AU or more. 

4-3. Binary mass ratios and very low-mass companions 

An estimate of the mass ratio for each binary candidate 
is given m Table 0; it has been obtained from the H band 
relative photometry, using th e 2M yr mass-luminosity re- 
lationship from Baraffc et al. ( 1998 ), which can be approx- 
imated by Mh = —3.25 log M + 2.19 for low-mass stars 
(M < IMq), and assuming that both stars are coeval 
and equally extincted. Because of the time dependency 
of the mass-luminosity relationship, these estimates are 
somewhat uncertain. We used the median age determined 
by H98 as a typical value for the whole cluster. Also, the 
extinction along the line of sight of the primary and the 
secondary are unknown and might be different (this effect 
should be rather small at 1.65 /zm, however), and infrared 
excesses can represent a significant part of the flux at this 
wavelength. In two cases (If A 144 and 157), the mass ratio 
estimated from A J is significantly smaller than from AH, 
indicating that at least one component shows a significant 
excess (IfA 143 and 158 secondaries actually show larger 
J — H excesses than their primaries by about 0.2 mag); 
in the other cases, both values are similar. Also, the H 
band relative photometry for IfA 184 is somewhat uncer- 
tain. However, we assume that these mass ratios are not 
systematically biased towards low or high values. 



We did not find any companion fainter than AH rs 
2 mag, although it would have been easily detected, pro- 
vided that the separation of the system is larger than 0'.'2 
(see Fig.|^). This flux ratio corresponds to a mass ratio of 
q = 0.25 at 2 Myr. For an average primary mass of 0.5M Q , 
this could point to the absence of very low-mass stars and 
brown dwarfs as secondaries. Alternatively, a statistical 
fluctuation cannot be excluded given the small number of 
detected binaries. 

In order to test the significance of this result, we per- 
formed Monte-Carlo simulations to compare the observed 
mass ratio distribution in IC 348 with that obtained as- 
suming that each companion has a mass (lower than its 
primary) taken at random from a given initial mass func- 
tion (IMF). Each simulated histogram is the average of 
1000 simulations, so that statistical uncertainties can be 
neglected; the exclusion of binaries with q < 0.1 does not 
modify the histograms by more than a few hundredths in 



each bin. First, we used the Kroupa et al. (1993) IMF, 
with ct\ = 1.3 and stellar masses in the range 0.08-1M©. 
The number of predicted binaries in the range q — 0- 



0.25 is about one. We also used the Reid & Gizis (1997a) 
IMF within the 5.2 pc solar neighbourhood, with and with- 
out brown dwarfs (i.e., with a minimum mass of 0.05 and 
0.075M Q respectively). In the first case, the mass function 
was chosed flat in the br own dw arf domain (ip(M) tx M°), 
following Reid & Gizis ( |l997b| ); we checked that the slope 
of the mass funtion does not modify significantly the re- 
sults. Without brown dwarf companions, we again predict 
about one companion in the first bin of Fig. |[ while this 
number is increased to almost three if we include substellar 
objects. In all three cases, the number of detected com- 
panions in the range q = 0.5-0.75 is about twice as large 
as that predicted by our simulations. In Fig. |], only bina- 
ries with Ma < 1M are plotted, because Reid & Gizis's 
IMF is only defined below IMq. It should be noted that 
the random pairing assumption together with the latter 
IMF leads to a mass-ratio distribution significantly differ- 



ent than that observed by Reid & Gizis (1997b ) . 

Despite of the small size of our small binary sample, 
the observations indicate that binaries in IC 348 are not 
preferentially equal mass systems, since no binary appears 
in the bin q = 0.75-1. According to DM91, the mass ratio 
distribution for solar-type field stars peaks near q = 0.2- 
0.3. On the other hand, in their studies of low-mass stars 
in the solar neighbourhood and the Hyades cluster, Reid 



& Gizis ( |l997a| , |l997b| ) concluded that M-dwarf binaries 
have a mass ratio distribution peaking at q w 1, which con- 
trasts with our findings fot IC348. It is unlikely that our 
observations have missed some equal mass systems; unless 
infrared excesses introduce a systematic bias against equal 
flux binaries, this deficit is real. 

The absence of small mass ratios (q < 0.25) in our 
observations is only marginally significant, given the sam- 
ple size. It is however more consistent with the models 
without brown dwarf companions. Reid & Gizis (1997b) 
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Fig. 4. Mass ratio distribution for binaries in IC 348 with 
Ma < ^Mq (thick histogram) compared to distributions 
simulated by taking at random the masses of t he com pan- 
ions from some IMFs. dotted: Kroup a et al. (1992) with 
«i = 1.3; dashed: Reid & Gizis ( 1997a ) in their 5.2 pc sam- 
ple, without brown dwarfs; dotted- dashed: Reid & Gizis 



(1997a) with brown dwarfs down to 0.05 M Q and with a 



flat mass distribution in the substellar domain. The latter 
model predicts a high number of binaries in the first bin, 
contrasting with the observations; it is very similar to the 
observed mass ratio distribution in the MS (DM91). The 
shaded histogram represents the subsample of low-mass 
stars (M < 0.3Af Q ) in IC 348. 



concluded from their Hyades study that the mass func- 
tion is flattening, maybe even decreasing, in the sub- 
stellar domain. However even a flat mass function prob- 
ably predicts too much companions with q < 0.25 in 
IC348. It rather seems that no brown dwarfs are com- 
panions of low-mass stars. This absence of brown dwarfs 
is also found in Table [j], where all companions have stellar 
masses, with the possible exception of the companion of 
HA192 (M B ~ 0.075M Q ). We caution again that these 
estimates are somewhat uncertain, due to age effects, in- 
frared excesses and, in some cases, unknown extinctions; 
spectra of the candidates should be obtained to determine 
their stellar/substellar status (IfA49, 82, and 184B have 
masses below 0.15 M Q from our results). 

If this absence of brown dwarfs is not due to statistical 
fluctuations or to systematic errors in the estimate of the 
mass ratio, this suggests that brown dwarfs cannot form 
in IC348 at separations larger than w 50 AU from stars. 
The results of Reid & Gizis ( 1997b| ) indicate that this limit 
is smaller than 5AU in the Hyades. This apparent lack 
of very low-mass companions may result from dynamical 
biasing during the early evolution of small subclusters: N- 



stars in a binary system and to the ejection of the lower 
mass components. However, in the solar neighbourhood, 
very low-mass secondaries (M < 0.1M©) can be found, 
at separations varying from 4 to 1800 AU (Reid & Gizis) 
1997a!) . 



5. Environmental conditions and binary 
formation 

In order to investigate evolutionary effects, we first com- 
pare the binary fraction we have determined in Sect.|^ 
with that of other SFRs (Sect. 5.1). The possible tempo- 
ral evolution of the binary fraction is discussed in Sect. |5. 2! 
before we argue that environmental conditions and binary 
frequency may be tightly linked. 



5.1. Comparison of binarity in IC 348 with other SFRs 

Given that the stars in IC 348 have a median age of 2 Myr, 
we can compare them directly with other SFRs like the 
Taurus- Auriga complex and the Orion Trapezium cluster, 
without introducing an evolutionary bias in terms of age. 
The binary excess observed in Taurus, Chamaaleon and 



Ophiuchus is of the order of a factor of 1.6 (Duchene 1998) 



A similar excess, if existing in IC 348, would yield a binary 
fraction of 37 ±5 % in our surveyed separation range. This 
is different from our observed binary fraction at a 2.5 a 
level, and can be excluded with a high confidence level 
(>98 likely harbor different binary fractions. From this we 
conclude that not all SFRs have a unique binary fraction, 
several of them exhibiting excesses (Taurus, Ophiuchus, 
Chamaeleon) and others showing binary fractions similar 
to that of the MS (the Trapezium cluster, IC 348, as well as 
other clusters in Orion: NGC 2024, 2068 and 2071, Ghez 



body simulations of Sterzik & Durisen (1998) show that, in 
most cases (~ 90%), the dynamical evolution of small- N 
systems result in the association of the two more massive 



et al. |1997| ). Consequently, we are led to the conjecture 
that the age, which is the same on average for the above- 
mentioned SFRs, is not the only parameter governing the 
binary fraction in a young cluster or a T-association. 

The stars we observed in IC 348 represent a rather 
large age spread, from a few 10 5 to about 10 7 years. This 
allows a comparison of the binary fraction with stellar age, 
which is shown in Table 0. We selected all stars in our 
samples with known age, which excludes the stars only 
detected in the Hcv survey, as well as two other members 
which lack I photometry in H98's study (If A 83 and 89). It 
appears that all three subsamples have similar binary frac- 
tions, although we are limited by the small sample sizes; 
a rank order test does not indicate any difference between 
the single and binary stars age distributions. This indi- 
cates that the binary fraction does not evolve significantly 
with time, at least on a timescale of a few Myrs. We also 
verified that the binaries in the three subsamples defined 
in Table ^ are roughly equally represented at all separa- 
tions. This means that we do not see any indication for 
an evolution of the binary separations over the timescale 
of the age range covered by these stars. 
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Table 5. Evolution of the binary fraction with primary 
age, as estimated by H98. There are a few stars in our 
sample with no age estimate (If A 83 and 89, and the Ha 
stars). 

t < 10 6 10 6 < t < 10 6 - 5 t > 10 6 ' 5 



observed 
binaries 



20 

5 



22 
3 



16 
4 



As an alternative to an evolutionary process, it has 
been proposed that the binary excess observed in vari- 
ous SFRs was the result of an observational bias: since 
the mass-luminosity relationship is shallower for younger 
stars, it is easier to detect companions around PMS stars 
than in the MS population (Zinnecker, priv. com.). The 
fact that we know several SFRs, now including IC348, 
with no binary excess indicates that this bias is not re- 
sponsible for the observed overabundance of PMS binaries 
in some regions. 

The binary fraction that seems to differ between SFRs 
relates to binaries which cover only a limited separation 
range. If the separation distribution was different from 
one SFR to another, the overall binary fractions could 
still be the same for all of them: an excess observed in a 
given separation range could be balanced by a deficiency 
of binaries with shorter or longer periods. There is cur- 
rently no such evidence, except perhaps for the study of 
the ROSAT population in Upper Scorpius by Brandner & 



Kohler (1998). Several arguments indicate that the orbital 
period distribution does not vary significantly between 
PMS and MS binaries: the number of spectroscopic bina- 
ries in Taurus is at least as large as that of the MS (Math- 
ieu 1994), lunar occultation surveys in this SFR have 
shown that the binary excess was present down to ~1 AU 
(Richichi et al. 1994; Simon et al. 



1995, and Pleiades bi- 



naries have a similar period distribution as dwarfs (Mer- 
milliod et al. 



1992; Bouvier et al. 1997) 



5. 2. Binary fraction and environmental conditions 



Considering IC348, the Trapezium (Prosser et al. 1994 



Petre et al. 1998) and Pleiades (Bouvier et al. 1997) clus- 
ters, and the solar neighbourhood stars (DM91), we have 
four samples with no binary excess at different evolution- 
ary stages (PMS, ZAMS and MS). We thus conclude that 
the binary fraction does not evolve with time between 
these stages. Any evolution of the binary frequency would 
have to occur within the first Myr after the formation pro- 
cess. Furthermore, the differing binary fractions between 
the various SFRs of the same age have to be explained, 
and a global time effect cannot be responsible for this. One 
common property of all the clusters without binary excess 
is that they are all rather dense: IC 348 has about 500 
stars pc~ 3 , and the Trapezium is about 10 times denser. 
The older Pleiades cluster, which is still dense nowadays, 



was probably even denser when younger, perhaps simi- 
lar to the Trapezium. On the other hand, the SFRs with 
high binary fractions (Taurus, Ophiuchus, Chamaeleon) 
are rather loose, with no more than a few stars pc~ 3 in 
the Taurus aggregates. This seems to indicate that a link 
exists between the binary fraction and the cluster density. 

Several physical processes could be the reason behind 
such a link. The impact of the average cluster density 
on the binary fraction could be direct; for instance, in 
dense clusters, the number of gravitational encounters is 
high and the binaries could be massively disrupted in such 
clusters over short timescales. From A-body simulations, 



Kroupa (1995a, 1995b) has shown that in clusters as dense 
as the Trapezium cluster the binary fraction could de- 
crease from 100% to about 50% in less than 1 Myr. Then, 
a model where all SFRs form with a high binary fraction 
(i.e., close to 100%), and where gravitational interactions 
between multiple systems are responsible for the decrease 
of the number of binaries, would be in qualitative agree- 
ment with the observational results: in all dense clusters, 
the binary fraction would have already decreased down to 
the MS level even for the youngest clusters in which bi- 
nary fractions have been measured so far, while it would 
have remained high in loose PMS associations. 

Alternatively, it is possible that the density is not the 
main parameter governing the binary fraction, but that 
another physical parameter, during or even before the 
star formation process, drives the subsequent evolution of 
the cluster, including its stellar density and binary frac- 
tion simultaneously. Durisen & Sterzik (1994) have shown 
that the current models of fragmentation and disk insta- 
bility may imply higher binary fractions when the initial 
cloud temperature is lower. In general, cold giant molecu- 
lar clouds may not be very efficient in forming stars, if the 
output of their fragmentation is small aggregates, with low 
densities (like Taurus). Regions creating high-mass stars, 
on the other hand, have rather high cloud temperatures; 
they usually from dense clusters, such as the Trapezium 
cluster. The link we find between binary frequency and 
cluster density could then be an intrinsic output of the 
fragmentation process. Other characteristics of the cloud 
before star formation occurs could as well be responsi- 
ble for the observed linked between density and binary 
fraction. For instance, the nature of the pre-collapse equi- 
librium in the parent cloud may influence the mass, size 
and angular momentum of the fragmented cores, leading 
to differing binary fraction and cluster density. 

Environmental conditions at the time of star formation 
thus could have an impact on the resulting binary popu- 
lation (i.e., the total number of multiple systems). At the 
present time, it is not possible to distinguish between a 
very rapid temporal evolution of the binary fraction, dur- 
ing the first Myr, or an intrinsic dependence of the binary 
fraction on these conditions. Observations of even younger 
populations in embedded clusters, as well as determination 
of accurate orbital period and mass ratio distributions are 
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needed to go further into the history of binary formation 
and evolution. 

6. Conclusions 

From a high-angular resolution study of IC 348 low-mass 
members wc find that the binary frequency in this very 
young cluster (~ 2 Myr) is similar to that of the Pleiades 
(~ 100 Myr) and of low-mass field dwarfs (~ lGyr). We 
therefore conclude that the binary frequency does not sig- 
nificantly evolve over time on a timescale of several 100 
Myr. Instead, it appears that the binary frequency among 
low-mass stars is already established at very young ages, 
i.e., within ~lMyr after the formation process. 

In particular, a long-term evolutionary effect cannot 
be responsible for the differing binary fractions found in 
different SFRs: on the one hand, Taurus and Ophiuchus 
exhibit binary excesses, on the other, the Trapezium and 
IC348 clusters don't. Yet, all these regions have similar 
ages of ~l-2 Myr. Furthermore, in the IC348 sample, we 
do not find evidence for an evolution of the binary fraction 
or the orbital period distribution within the age spread of 
the cluster of about 10 Myr. A time evolution of the bi- 
nary frequency, if any, is thus constrained to occur within 
the first 1 Myr of stellar evolution. After this time, intrin- 
sic differences exist between SFRs regarding their binary 
content. 

In spite of the large dynamic range of our images, no 
brown dwarf companion is found in IC 348 binaries (to 
the possible exception of KA192B). The mass ratio dis- 
tribution we find is consistent with the absence of brown 
dwarf companions to low-mass members of this cluster; 
a similar conclusion was drawn for the Hyades cluster. 
Also, the mass-ratio distribution is not peaked towards 
q = 1, in possible contradiction with what has been pro- 
posed for the solar neighbourhood M-dwarfs population. 
Further studies with larger telescopes will allow a better 
determination of this distribution. 

Comparing the results obtained on IC 348 to similar 
studies in other clusters, it appears that the binary frac- 
tion may be inversely correlated with the average clus- 
ter density, with dense clusters showing low binary frac- 
tions (similar to field dwarfs), as opposed to the loose T- 
associations like the Taurus- Auriga and Chamacleon com- 
plexes where the binary fraction is larger. On the basis 
of this qualitative trend, at least two scenarios may ex- 
plain the observed differences in binary fractions: either 
the formation mechanism always leads to an initially high 
binary fraction (of the order of 100%) and frequent gravi- 
tational encounters in dense clusters disrupt binaries on a 
timescale of 1 Myr or less, or specific initial conditions in 
the parent molecular cloud, such as the gas temperature, 
metallicity, angular momentum , etc., lead to different 
output of the star formation process and govern simulta- 
neously the binary frequency and the cluster density. High 
resolution studies of embedded clusters even younger than 



those investigated so far are still needed to settle these is- 
sues. 
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